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The description of the proximity effect in superconducting/ferromagnetic heterostructures requires 
to use spin-dependent boundary conditions. Such boundary conditions must take into account the 
spin dependence of the phase shifts acquired by electrons upon scattering on the boundaries of fer- 
romagnets. The present article shows that this property can strongly affect the critical temperature 
and the energy dependence of the density of states of diffusive heterostructures. These effects should 
allow a better caracterisation of diffusive superconductor/ferromagnet interfaces. 

PACS numbers: 73.23.-b, 74.20.-z, 74.50.+r 
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I. INTRODUCTION 

When a ferromagnetic metal (F) with uniform mag- 
netization is connected to a BCS superconductor (S), 
the singlet electronic correlations characteristic of the 
S phase can propagate into F because electrons and 
holes with opposite spins and excitation energies are cou- 
pled coherently by Andreev reflections occurring at the 
S/F interface. Remarkably, the ferromagnetic exchange 
field induces an energy shift between the coupled elec- 
trons and holes, which leads to spatial oscillations of the 
superconducting order parameter in F— This effect 
has been observed experimentally through oscillations 
of the density of states (DOS) in F with the thickness 
of F— , or oscillations of the critical current Jo through 
S/F/S structures^ 5 -' 6 -' 7 -, with the thickness of F or the 
temperature. The oscillations of Io have allowed to ob- 
tain 7r-junctions£, i.e. Josephson junctions with I < 0, 
which could be useful in the field of superconducting 
circuits 9,10 . A reentrant behavior of the superconducting 
critical temperature of S/F bilayers with the thickness 
of F has also been observed-^. At last, some F/S/F 
trilayers have shown a lower critical temperature for an 
antiparallel alignment of the magnetizations in the two F 
layers as compared with the parallel alignment^, which 
should offer the possibility of realizing a superconducting 
spin-switc h 13 : 14 . 

For a theoretical understanding of the behavior of 
S/F hybrid circuits, a proper description of the inter- 
faces between the different materials is crucial. For a 
long time, the only boundary conditions available in 
the diffusive case were spin-independent boundary con- 
ditions derived for S'/normal metal interfaces^ 5 . Re- 
cently, spin-dependent boundary conditions have been in- 
troduced for describing hybrid diffusive circuits combin- 
ing BCS superconductors, normal metals and ferromag- 
netic insulators^. These boundary conditions take into 
account the spin-polarization of the electronic transmis- 
sion probabilities through the interface considered, but 
also the spin-dependence of the phase shifts acquired by 
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FIG. f: a. Diffusive F/S/F trilayer consisting of a BCS 
superconductor S with thickness ds placed between two fer- 
romagnetic electrodes Fi and 7*2 with thickness d_F- In this 
picture, the directions of the magnetic polarizations in Fi and 
F-2. are parallel [antiparallel], which corresponds to the con- 
figuration C = P [AP]. b. S/F bilayer consisting of a BCS 
superconductor S with thickness ds/2 contacted to a ferro- 
magnetic electrode F with thickness dp. 



electrons upon transmission or reflection by the interface. 
The first property generates widely known magnetoresis- 
tance effects 1 ^. The second property is less commonly 
taken into account. However, the Spin-Dependence of In- 
terfacial Phase Shifts (SDIPS) can modify the behavior 
of many different types of mesoscopic circuits with ferro- 
magnetic elements, like those including a diffusive normal 
metal island 1 ^, a resonant syste m 19 ' 20 , a Coulomb block- 
ade syste m 20 ' 21 ' 22 , or a Luttinger liquid^ 3 -. It has also 
been shown that the SDIPS has physical consequences 
in S/F hybrid system s 16 : 24 ' 25 : 26 . One can note that, in 
some references, the SDIPS is called "spin-mixing angle" 
or "spin-rotation angle" (see e.g. Refs. [24ll26h . In the 
diffusive S/F case, the spin-dependent boundary condi- 
tions of Ref. have been applied to different circuit 
geometrie s 27 ' 28 ' 29 ' 30 ' 31 but the only com pari son to exper- 
imental data has been performed in Ref. [29|. The authors 
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of this Ref. have generalized the boundary conditions of 
Ref. [ID to the case of metallic S/F interfaces with a su- 
perconducting proximity effect in F. They have showed 
that the SDIPS can induce a shift in the oscillations of 
the critical current of a S/F/S Josephson junction or of 
the DOS of a S/F bilayer with the thickness of F. Signa- 
tures of this effect have been identified in the Nb/PdNi 
hybrid structures of Refs. Nevertheless, the problem 
of characterizing the SDIPS of diffusive S/F interfaces 
has raised little attention so far, in spite of the numerous 
experiments performed. 

A good characterization of the properties of diffusive 
S/F interfaces would be necessary for a better control 
of the superconducting proximity effect in diffusive het- 
erostructures. The present article presents other conse- 
quences of the SDIPS than that studied in Ref. 1 2 9. which 
could be useful in this context. In particular, the SDIPS 
can generate an effective magnetic field in a diffusive S in 
contact with a diffusive F, like found for a ballistic S in 
contact with a ferromagnetic insulator 24 . This effective 
field can be detected, in particular, through the DOS of 
the diffusive F layer, with a visibility which depends on 
the thickness of F. A strong modification of the varia- 
tions of the critical temperature of diffusive S/ F struc- 
tures with the thickness of F is also found. These effects 
should allow to characterize the SDIPS of diffusive S/F 
interfaces through DOS and critical temperature mea- 
surements, by using the heterostructures currently fabri- 
cated. The calculations reported in this paper are also 
appropriate to the case of a diffusive S layer contacted 
to a ferromagnetic insulator (FI). 

This paper is organized as follows: Section II presents 
the initial set of equations used to describe the het- 
erostructures considered. The case of F/S/F trilayers 
is mainly addressed, but the case of S/F (or S/FI) bi- 
layers follows straightforwardly. Section III specializes 
to the case of a weak proximity effect in F and a super- 
conducting layer with a relatively low thickness ds < £s, 
with £s the superconducting coherence length in S. The 
spatial evolution of the electronic correlations in the S 
and F layers is studied in Section III. A. The energy- 
dependent DOS of S/F heterostructures is calculated in 
Section III.B. Section III.C considers briefly the limit of 
S/FI bilayers. Section III.D discusses SDIPS-induced 
effective field effects in other types of systems. Section 
III.E compares the present work to other DOS calcu- 
lations for data interpretation in S/F heterostructures. 
Critical temperatures of S/F circuits are calculated and 
discussed in Section III.F. Conclusions are presented in 
Section IV. Throughout the paper, I consider conven- 
tional BCS superconductors with a s-wave symmetry. 



II. INITIAL DESCRIPTION OF THE PROBLEM 

This article mainly considers a diffusive F / S/F tri- 
layer consisting of a BCS superconductor S for — ds/2 < 
x < ds/2, and ferromagnetic electrodes F\ for x <E 



{-d s /2 - d F , -d s /2} and F 2 for x G {d s /2, d s /2 + d F } 
(see Figure I. a). The magnetic polarization of the two 
ferromagnets can be parallel (configuration C = P) or 
anti-parallel (configuration C = AP), but the modu- 
lus \E ex \ of the ferromagnetic exchange field is assumed 
to be the same in F\ and F%. Throughout the struc- 
ture, the normal quasiparticle excitations and the su- 
perconducting condensate of pairs can be characterized 
with Usadel normal and anomalous Green's functions 
G„, CT = sgn(a; n ) cos(0„ jO -) and F n>a = sin(0„ !<7 ), with 
O n ,a-(x) the superconducting pairing angle, which de- 
pends on the spin direction a G {T,l}, the Matsubara 
frequency uj n {T] = (2n + VjirksT, and the coordinate x 
(see e.g. Ref. 1321). The Usadel equation describing the 
spatial evolution of 9 n ^ writes 



hDs d 2 e n ^ 

2 dx 2 
in S and 

hD F d 2 fl„ |g 

2 dx 2 



\w n \ sin(0„ )(7 ) - A(x) cos(0 n>a ) (1) 



(|w„ + iE ex asgn(uj n )) sm(9 n ^) (2) 



in Fi and F 2 , with Dp the diffusion constant of the fer- 
romagnets and Ds the diffusion constant of S. The self- 
consistent superconducting gap A(x) occurring in JT]) can 
be expressed as 

A(z) = ^!^ £ sin(0„ >CT ) (3) 
^e{T,l} 

w„(T)£{-Q D ,n D } 

with p> the Debye frequency of S, A -1 = 
27rk B T c BCS J2 Un(T sc S)e{0tQD} u- 1 the BCS coupling 

constant and T BGS the bulk transition temperature of 
S. I assume A = in F\ and F 2 . The above equations 
must be supplemented with boundary conditions describ- 
ing the interfaces between the different materials. First, 
one can use 



d0„ 



dx 



o 



(4) 



x=±(d s /2+d F ) 



for the external sides of the structure. Secondly, the 
boundary conditions at the S/F interfaces can be cal- 
culated by assuming that the interface potential locally 
dominates the Hamiltonian, i.e. at a short distance it 
causes only ordinary scattering (with no particle-hole 
mixing) (see e.g. Ref. 33)- This ordinary scattering 
can be described with transmission and reflection am- 
plitudes tn,a and rn,a ^ for electrons coming from the 
S(F) side of the interface in channel n with a spin 

direction a. The phases of tn^a and rn,a can be 
spin-dependent due to the exchange field E ex in Fi(2) 
and a possible spin-dependence of the barrier poten- 
tial between S and ^1(2) ■ Boundary conditions tak- 
ing into account this so-called Spin-Dependence of In- 
terfacial Phase Shifts (SDIPS) have been derived for 
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|i^ iT | 2 , |tfj 2 < 1 and a weakly polarized Fi£i22. When 
there is no SDIPS, the boundary conditions involve the 
tunnel conductance Gt = Gq J2 n Tn and the magne- 
toconductance G M n, = Gq ^„(|^ >t | 2 - I^J 2 ), witn 
t (J.) the majority(minority) spin direction in the F elec- 
trode considered, Gq = e 2 /h, and T n = |i.^ T | 2 + l^jj 2 - 
In the case of a finite SDIPS, one must also use the 

F(S) _ on r f F(S) A ,S(F) 



conductances G^ {J ' — 2GqJ2 u (p 



G 



4[<^VT n ]), 

G Q J2 n Tn(pn iS) + 

T ? ? (F) )/4, with p™ = Imfr^r^*] and r™ = lm[t^ L *} 
for to £ {S, F}. In the following, I will focus on the ef- 



F(S) 

G x to be negligible, like found with a simple barrier 
model in the limit T n < 1 and E ex <C Ep^. In this 
case, one finds that the boundary conditions for the S/F 
interface located at x = s 
write 



fccts of and Ga, and I will assume Gmr, G^ S ^ and 



(-l) J d fl /2, with j G {1,2}, 



0# 



(-iy 7T sin[< 



9 s 1 

n,cri 



(5) 



and 



00* 



<9:r 



7 



E 



(6) 

where the indices 5 and F indicate whether 9 nM 
and its derivative are taken at the S or F side 
of the interface. These equations involve the re- 

F( S) 

duced conductances jt = Gt£,f /Aap and 7a = 

F( S) 

^f/Agf, the barrier asymmetry coefficient 7 = 
£,s< j f/£,f&s, the superconducting coherence lengthscale 
£s = (hDs/ZABcs) 1 ^ 2 > the magnetic coherence length- 
scale £f = (?iDf/ \E ex \y/ 2 , the gap Abcs for a bulk 5, 
the normal state conductivity o~prs) of the F (<S) mate- 
rial, and the junction area A. The coefficient j takes 
into account the direction of the ferromagnetic polar- 
ization of electrode Fj in configuration C G {P, ^4P}. 
One can use the convention e„'£ = (— iy asgii(uj n ) and 
e nCT ,J = (JS g rL ( a -'ri)' m which the factor sgn(w n ) arising 
from the definition of Q n ^ a and the terms (— 1) J and a 
arising from the boundary conditions have been included 
for compactness of the expressions. Note that in the pres- 

F( S) 

ence of a finite SDIPS i.e. 7a ^ 0, the right hand side 
of equation © is not zero contrarily to what found in 
the spin-degenerate case^. In the general case, 7^ and 
7^ are different (see e.g. Appendix A). This implies that, 
with the present approximations, each interface is char- 
acterized by three parameters: 77-, 7^ and 7?. For the 
sake of simplicity, symmetric F/S/F trilayers are con- 
sidered, so that 7t, 7^ and 7^ are the same for the two 
S/F interfaces. 

Before working out the above system of equations, it is 
interesting to note that the angle Q n ^ calculated in the 



parallel configuration C — P for x > also corresponds 
to the angle 9 n , a expected for a S/F bilayer consisting of 
a superconductor S for < x < ds/2, and a ferromag- 
netic electrode F for x £ {ds/2, ds/2 + dp} (Figure l.b). 
In practice, using a F/S/F geometry can allow one to 
obtain more information on spin effects, as shown below. 



III. CASE OF A THIN SUPERCONDUCTOR 
AND A WEAK PROXIMITY EFFECT IN F 

A. Spatial variations of the pairing angle 

I will assume that the amplitude of the superconduct- 
ing correlations in Fu 2 ) is weak, i.e. \6 n ,<r\ ^ 1 f° r 
x G {-ds/2 - d F , -d s /2} and x G {d s /2,d s /2 + d F } 
(hypothesis 1) so that one can develop the Usadel equa- 
tion @ at first order in n<cr . This leads to 



dx 2 



'k c ^ 
1f~ 



= 



(7) 



in the ferromagnet Fj, with j G {1,2}, k^'i = {2[ir)^ a + 

\uJ n /E ex \)]) 1/2 and fftj, = (-lpef^. Combining Eqs. 
(0| and (7}, one finds in Fj 



cosh [x - (-iy (d F + 4f )] 



cosh ( dp- 



(8) 

This result together with boundary condition © leads 
to 



^n,(j{ x j) 



7TSin(^ (xj)) 



7T cos(0| i<y (ay)) + + B% 



(9) 



with B^'l = k^'i- tanh[dirfc^'^./^ir]. This allows to rewrite 
the boundary condition ([6]) in closed form with respect 



to 



80 s 



dx 



(-1)^4^ sin^^.)] (10) 



for the S/F interface located at x — Xj, with 



J.n,<T 



7 7 t cos(0£ CT ) + + I-!* if n % 

In the following, I will assume \E ex \ 3> Abcs like m most 
experiments, so that k^ a = 1 + ir\^ a - I will also assume 
ds/^S 5~ 1, so that one can use, for the C configuration 
and -ds/2 < x < ds/2, 



71, a n,(T\ 



<>/&)" /3n, ff (*/£s) 2 (12) 



4 



with 



<C 1 {hypothesis 2). Note that al- 
though experiments are often performed in the limit of 
thick superconducting layers ds > £s, assuming d$ < £s 
is not unrealistic since one can obtain diffusive supercon- 
ducting layers with a thickness ds ~ £s (see e.g. Ref . l34l) . 
Furthermore, using relatively low values of ds is more 
favorable for obtaining efficient superconducting spin- 
switches^. Hypothesis 2 allows one to develop sin(#„ j(T ) 

and cos(0 n)(7 ) at first order with respect to 9^-^ ff in S. 
Accordingly, 1 will neglect the space-dependence of A{x) 
and assign to it the value A c in configuration C. Note 
that I do not make any assumption on the value of the 
angle 0% a , which is not necessarily close to the bulk BCS 
value. The coefficient Cj na is transformed into its conju- 
gate when the magnetic polarization of electrode Fj is re- 
versed. Therefore, I will note £f n a — C%*£„ — C„ and 
£i n a = ifinf • ^ ne & bove assumptions lead to a^ a = 0, 



P 

71,(7 



4£° 



4^ + 003(^^)51^ 



(13) 



4sm(9i a ) ~ cos(9i a )6 s (3. 



AP 



and 



B AP 



4 + 2Re[£-]cos(^ )CT )5 s 
(8Re[££] +4|£-| 2 cos(^ J^) sin^ J 



(14) 



S^ + eRe^cos^J^ + |£-|^cos 2 (^J5| 

(15) 

with <5s = ds/£,s- On the other hand, from ([1]), one finds 

flC A c cos(^,J~|^|sin(^ CT ) 

= (16) 



The comparison between equations (fill)) . (fT5)) and (pi 
allows one to find 9^ a as a function of A c . Then, one has 
to calculate A c by using the self-consistency relation |[3J • 
I will study below the DOS and the critical temperature 
following from these Eqs., in a limit which leads to simple 
analytical expressions. 



B. Low-temperature density of states of S/F 
heterostructures 



The DOS of the ferromagnets F\ and F% of Figure l.a 
can be probed at x = i{dp + dg/2) by performing tun- 
nelling spectroscopy through an insulating layer. So far, 
this quantity has been less measured 3 ^ 36 ' 37 ! 38 than criti- 
cal temperatures or supercurrents. However, this way of 
probing the superconducting proximity effect is very in- 
teresting because it allows one to obtain spectroscopic 
information. It has been shown that the zero-energy 
DOS of a F layer in contact with a superconductor os- 
cillates with the thickness of F. For certain thicknesses, 



this zero-energy DOS can even become higher than its 
normal state value Nq— i 40 i 41 i 42 , as shown experimen- 
tally in Ref. H Remarkably, the SDIPS can shift these 
oscillations 2 -. Although the energy dependence of the 
DOS of diffusive S/F structures has raised some theoret- 
ical and experimental interest, the effect of the SDIPS on 
this energy dependence has not been investigated so far. 

For calculating analytically the low-temperature 
DOS of the structure of Fig. l.a., one can as- 
sume cos(6»f (T (a; i ))/ |( j T cos{9^ tJ {x j )) 



■ F P 2 



B n'Z){lT + ij^Vnla )l < 1 {hypothesis 3), which leads 
to L a n — 7 (jt + *70 cr sgn(o;„)^ . This hypothesis is e.g. 

valid for dp > £p and any value of Jj, if 7t is rel- 
atively small (see e.g. Fig. [2]). I will also assume 
that the lowest order terms in 5s prevail in the numera- 
tors and denominators of expressions (| 1 3[) and (|15p. i.e. 
/3£ ff ~ Oin^)^ 1 and ^ ~ Re[££] flin^)^ 1 
{hypothesis 4-)- Taking into account hypothesis 3 and 



-4 



are 



7 ~ 1, hypothesis 4 is valid provided 7t and 
relatively small compared to 1. Importantly, hypotheses 
3 and 4 are less restrictive regarding the value of 7 J . Ac- 



cordingly, I will often use 7t, 7^ <C 
this paper. Hypotheses 3 and 4 lead to 



= arctan 



A c 



Mr. 



+ n c 

1 n, tr 



in the Figs, of 



(17) 



with 



and 



= 2A BC s7(7t + i7^sgn(w„))<S5 1 (18) 



AP 



2&bcs11t5s 



(19) 



Equations (fT7|) and (fT8|) show that an effective magnetic 
field H e f f appears in the S layer in the P configuration, 
due to 7^ ^ 0. From Eqs. {JXTP and {HJ), H eff can be 
expressed as 



gHBH e ff 



d s 3a s A 



qS 

2E TH^ (20) 



Here, v§. and if denote the Fermi velocity and mean 
free path in S, and Gs — asA/ds and E§, H = HDs/dg 
denote the normal state conductance and the Thouless 
energy of the S layer of Fig. l.a. From Equations (JTTJ) 
and (TT9"]) . the effective field effect disappears in the AP 
configuration because the two contacts are assumed to 
be symmetric, and therefore, their contributions to -He// 
compensate each other in the AP case. Note that, in 
principle, the 7^ term can induce an effective field ana- 
logue to H e ff in the F layer, but this effect is not rele- 
vant in the regime studied in this paper (see Appendix 
B). The effects of H e ff on the DOS of the structure will 
be investigated in next paragraphs. In order to calculate 
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y=1,y T =0.12, d s /^ s =1, case P 




Y ; = (A P = 0.72V S ) 
1,^ = (A^O.72^) 
2.1, v s = 0.1 2 (A P = 0.70^) 
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FIG. 2: Zero energy density of states N (e = 0) at x = 
<^f + ds/2 as a function of c?f, for the F/S/F structure of 
Figure l.a in the P configuration, with 7 = 1, jr = 0.12, 
and ds/^s ~ 1- The different curves correspond to different 
SDIPS parameters: j F = 7^ = (black dotted curve), {7^ = 
-2.1, 7| = 0} (black full curve), and {7^ = -2.1, 7 f = 0.12} 
(red full curve). The low-temperature self-consistent gap A p 
found for these three different cases is indicated at the bottom 
right of the Figure. Using a strong value for -y F allows one 
to change significantly the phase of the oscillations of N F (0) 
with d,F (effect similar to Ref. [29j). One can check that the 
hypotheses 1 to 4 are valid for the parameters used in this 
Figure. Note that for the different cases considered here, the 
critical temperature T P of the structure (see Sec. III.F) is 
such that 0.615T C SCS < T P < 0.635T FCS . The four red 
points correspond to the red full curves shown in Fig. [3] 

one has to combine the self-consistency relation ([3]) 
with Eq. (|17p . which gives, at low temperatures, 



Re[log( 



(A c r 



A 



BCS 



)] = (21) 



This equation can be solved numerically. The resulting 
A c is independent from the values of n and a used in 
Eq. (|2ip. Then, the value of the pairing angle 6 n . a in 
the ferromagnets can be found by using Eqs. ©, (O, 
dHJ) and (17]). Note that for 7 < 1, the above Eqs. 
are in agreement with formula (5) of Ref. [29|, obtained 
with rigid boundary conditions, i.e. 9 n ^ a equal to its bulk 
value at the S side. The energy dependence of rha can 
be found by performing the analytic continuation oj n — 
—is + T and sgn(oj„) = 1 in the above equations. The 
rate F = 0.05 is used to account for inelastic processes^. 
At last, the density of states N a (x,e) at position x for 
the spin direction a £ {|, j} can be calculated by using 
N a (x,e) = (iVo/2)Re[cos[0 n , (T (x)]], where N /2 is the 
normal density of states per spin direction. 

In the following, I will mainly focus on N F (e) = 
Scre{t i}^ cr ( x = dF,s)- Figure shows the variations 
of the zero energy density of states N F (e = 0) as a func- 
tion of dp, for interface parameters 7 t = 7^ = 0.12 and 

7^ = —2.1. Importantly, the value 7 t = 0.12 seems real- 
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FIG. 3: Density of states N F (e) at x = dF+ds/2 as afunction 
of the energy e, for the F/S/F structure of Figure l.a with 
7 = 1, 7 t = 0.12, ds/is = 1, 7 J = -2.1 and 7 £ = 0.12. 
The different panels correspond to different values of o!f/£f- 
The red full curves correspond to the P configuration and the 
blue dashed curves to the AP configuration. For finite values 
of some characteristic "double-gap" structures appear in 
N F (e) in the P case, due to a SDIPS-induced effective field 
appearing in S. The visibility of this effect strongly depends 
on the thickness dF of the F layers. For the different cases 
considered here, one can check that hypotheses 1 to 4 are 
valid and that the critical temperature T% of the structure 
(see Sec. III.F) is such that 0.62T C SCS < T c c < 0.67T C SCS . 
Note that for the above parameters and 7^ = 0, the curves 
obtained in the P and AP case would be identical and very 
close to the AP curves shown here. 



istic, at least for the weakly polarized Nb/Pdo.gNio.i bi- 
layers used in Ref. [3|, for which one finds jt ~ 0.15 (see 
Ref. [lil ). In addition, a simple barrier model suggests 



that the situation 



G 



3> Gt, with G F < can happen 
(see Appendix A). The value 7^ = —2.1 used in Fig. [3] 
thus seems possible. One can see that j F can change 
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significantly the phase of the oscillations of N F (e = 0) 
with dp (this effect has already been studied in Ref. [23 
in the case of rigid boundary conditions but I recall it 
here for the sake of completeness). In Fig. [21 using for 
7^ a strong negative value allows to get N F (0) < No for 
the lowest values of d,F , like often found in experiments. 
Note that from Fig. [5[ 7^ can also shift the oscillations 
of N F (e = 0) with dp- Here, this effect is much weaker 
than that of j F , but one has to keep in mind that the 
limit 7^ <C J F is considered. 

Equation (|20| shows that a measurement of H e ff 
should allow to determine the conductance of a diffu- 
sive S/F interface. In this context, studying the energy 
dependence of N F (e) is very interesting, because it can 
allow to see clear signatures of FI e ff, as shown below. 
Figure [3] shows the energy dependence of N F (e) in the 
P and AP configurations, for a finite value of 7^ and 
different values of cLf- For C — P, N F (e) shows some 
"double-gap" structures which disappear if the device is 
switched to the AP configuration. These double struc- 
tures are an indirect manifestation of the effective mag- 
netic field H e ff occurring in S in the P configuration, 
due to 7^ / 0. Although i? e // is localized in the S 
layer, the double-gap structure that this field produces 
in the DOS of S is transmitted to the DOS of F due to 
Andreev reflections occurring at the S/F interfaces, as 
shown by Eq. ©. Interestingly, Rowell and McMillan 
have already observed that an internal property of a S 
layer can be seen through the superconducting proximity 
effect occurring in a nearby normal layer. More precisely 
these authors have found that the DOS of an Ag layer 
can reveal the phonon spectrum of an adjacent supercon- 
ducting Pb layer~. Remarkably, the visibility of H e t t in 
N F (e) is modulated by quantum interferences occurring 
in F. Indeed, H e ff is more visible for certain values of 
d,F (e.g. (If/£f = 1-0 or 1.2 in Fig. [3j) than others (e.g. 
df/^F = 2.1 in Fig. [3]), due to the dF-dependence of 
Eq.©. 

It is useful to note that the SDIPS-induced effective 
field H e f f should also occur in the S/ F bilayer of Figure 
Lb. In this case, the Thouless energy and normal state 
conductance of the S layer correspond to Ej, H = iEj, H 
and Gs — 2Gs respectively, so that one finds gusHeff = 
Ej, H G^/Gs- Double gap structures strikingly similar to 
those shown in Figure [3] were indeed measured very re- 
cently by P. SanGiorgio et al., at the ferromagnetic side 
of diffusive Nb/Ni bilayers, in the absence of any exter- 
nal magnetic field 4 ^. Remarkably, the visibility of the 
observed double structures varies with dp, as predicted 
above. Note that in S/F bilayers, the field H e ff should 
also be observable directly at the S side by measuring 
N s (e) — J^o-eft 4} ^<r( x ~ 0i e )- However, for param- 
eters comparable to those of Figure [3J this should not 
enhance the resolution on H e ff (see figure [H left). For 
certain values of dp, H e ff is even more visible in N F (e) 
than in N s (e)(see Figure UJ. 

Before concluding this section, I would like to empha- 
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FIG. 4: Energy dependent density of states N s (e) measurable 
at the S side (x = 0) of the S/F structure of Figure l.b. The 
left panel corresponds to the case of a metallic F contact 
with parameters corresponding to that of Figs. [2] and |3j In 
this case, the DOS measured at the S side does not allow to 
resolve the SDIPS-induced effective field H e jj better than a 
DOS measurement at the F side of the structure, as can be 
seen from a comparison with Fig. [3] For certain values of (If, 
Heff is even more visible in N F (e) than N s (e). The right 
panel corresponds to the case in which F is not a metal but 
an insulating ferromagnet, i.e. 7t = . In this case, one 
can use the reduced SDIPS parameter = Gf^s/Aas- For 
7^ 0, N s (e) shows strong signatures of the effective field 
Heff induced by the FI layer in S. One can check that the 
hypotheses 1 to 4 are valid for the parameters used in this 
Figure. 

size that from Eqs. (fT?)) and ([r5|) , in the AP config- 
uration, the SDIPS-induced effective field H e ff disap- 
pears for the F/S/F structure considered in this paper 
because the two contacts are assumed to be symmetric 
and have thus opposite contributions to H e ff in the AP 
case. In the case of a dissymmetric structure, this should 
not be true anymore, but the SDIPS-induced effective 
field should nevertheless vary from the P to the AP case. 
This is one practical advantage of working with F/S/F 
trilayers instead of S/F bilayers. 

C. Low-temperature density of states of S/FI 
bilayers 

Twenty years ago, internal Zeeman fields were ob- 
served in superconducting Al layers contacted to different 
types of ferromagnetic insulators (FI) (see Refs. I49ll50l. 
I51H52T ). Using a ballistic S/FI bilayer model, Ref. [2J 
suggested that the observed internal fields could be in- 
duced by the SDIPS 24 . However, the inadequacy of this 
ballistic approach for modeling the actual experiments 
was pointed out in Ref. [52j. Most of the experiments 
on AX/FI interfaces were interpreted by their authors 
in terms of a diffusive approach with no SDIPS, and an 
internal Zeeman field added arbitrarily in the Al layer 
(see Refs. [UliiiS). The calculations of Section III.B. 
provide a microscopic justification for the use of such an 
internal field in the diffusive model. Indeed, using "fx — 
in the above calculations allows one to address the case 
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of diffusive S/FI bilayers. One finds that the SDIPS- 
induced effective field H e ff of Eq. (J^DJ) can occur in a 
thin diffusive S layer contacted to a FI layer. This effec- 
tive field effect can be seen e.g. in the density of states 
N s (e) of the S layer at x = (sec Figure Q] right). Re- 
markably, it was found experimentally 52 that H e f / scales 
with dj 1 , in agreement with Eq. (f2TJ)^ . 

D. SDIPS-induced effective fields in other types of 
system 

Interestingly, the SDIPS can induce effective field ef- 
fects in other types of systems. First, the case of S/N/FI 
trilayers with a thickness djv of normal metal N has 
been studied theoreticall y 27 ! 28 . In this case, a con- 
ductance G^ similar to G^ can be introduced to take 
into account the SDIPS for electrons reflected by the 
FI layer. The N layer is subject to an effective field 
g^BH' e ff — EthG^ /Gn with Gn the conductance of N. 
The expression of H' e ff is analogue to that of H e f / (see 
Eq. 120)1 . up to a factor 2 which accounts for the symme- 
try of the F/S/F structure with respect to x = in the 
P configuration. Secondly, an effective field H e ff defined 

by gVBH'eff = (Kvf /2L)((fl + ip\ -<p l L - <Pr) ha s been 
predicted for a resonant single-channel ballistic wire with 
length L placed between two ferromagnetic contacts^, 
with the reflection phase of electrons with spin a 

incident from the wire onto the left(right) contact, and 
is the Fermi velocity in the wire. The expression of 
H e fj also shows strong similarities with that of Heff ( see 
Eq. [20j middle term). In practice, signatures of the field 
H e ff could be identified in a carbon nanotube contacted 
with two ferromagnetic contacts 2 - 2,56 . The fields H e ff, 
H e ff and H e ff have the same physical origin: the en- 
ergies of the states localized in the central conductor of 
the structure depend on spin due to the spin-dependent 
phase shifts acquired by electrons at the boundaries of 
this conductor. In all cases, the DOS of the central con- 
ductor reveals the existence of the SDIPS-induced effec- 
tive field only if it already presents a strong energy depen- 
dence near the Fermi energy in the absence of a SDIPS. 
In the F/S/F case, this energy dependence is provided 
by the existence of the superconducting gap in S*L. In 
the S/N/FI case, it is provided by the existence of a su- 
perconducting minigap in TV. At last, in the case of the 
ballistic wire, it is provided by the existence of resonant 
states in the wire. 



E. Comparison between the present work and 
other models for data interpretation in S/F 
heterostructures 

For characterizing the properties of S/ F interfaces, one 
has to interpret the experimental data showing the oscil- 
lations of the density of states N F (e) in F with the thick- 



ness dp of F (or the oscillations of the critical current Io 
of a S/F/S Josephson with dp)- However, if one uses a 
simple description with spin-degenerate boundary condi- 
tions, the amplitude and the phase of these signals are 
not independent, which makes the agreement with ex- 
perimental results impossible in most cases. The SDIPS 
concept can solve this problem since it produces a shift 

S(F) 

of the signals oscillations with respect to the = 

case. However, in many cases, the observed shifts were 
attributed to the existence of a magnetically dead layer 
(MDL) at the F side of the S/F interface (see e.g. Refs. 
HHUH). In other cases, the discrepancy between the 
theory and the data was solved by taking into account 
spin-scattering processes in the F layer (see e.g. Ref. @). 
In order to have a better insight into superconducting 
proximity effect experiments, one must stress the impor- 
tance of estimating experimentally the MDL thickness 
and the spin-scattering rate. In principle, spin-scattering 
rates can be estimated experimentally, as was done for 
instance for the CuNi alloy 6 ^ which is frequently used 
in proximity effect measurements (see e.g. [6. 58). An 
experimental determination of the MDL thickness has 
also been performed in a few structures used to measure 
T c or lySLSMZ&MS., but, so far, this parameter has not 
been used for a real quantitative analysis of the data. In 
some situations, a model combining the SDIPS with spin- 
scattering and/or a MDL may be necessary. In any case, 
it is important to point out that descriptions based on 
spin-degenerate boundary conditions are, in principle, in- 
complete since they do not account for the effective field 
effect described in section III.B. 

Before concluding this section, it is interesting to note 
that the effective field effect produced by G? in S or the 
phase shift of the spatial oscillations of the DOS provided 
by G F in F will remain qualitatively similar when the 
signs of G^ and G F are changed (not shown). The sign 

of G F in the experiments of Refs.dlll could be determined 
from a quantitative comparison between the theory and 
the data 2 ^ . Below, we present a study of the critical tem- 
perature of S/ F structures which can give information on 
the signs of G? and G F through qualitative signatures. 



F. Critical temperature of S/F heterostructures 

The critical temperature of S/F hybrid structures 
has already been the topic of many theoretical (see 
e.g. Refs. [HHl 4,35,6611671) and experimental (see e.g. 
Refs. [Illlil) studies, but the effects of the SDIPS on this 
quantity have raised little attention so far. I show below 
that the SDIPS can significantly modify the critical tem- 
peratures of S/F diffusive structures. Calculating the 
critical temperature Tjf of the structure of Figure I. a in 
configuration C requires to consider the limit in which 
superconducting correlations are weak in S as well as in 
F (hypotheses 1 and 2 are then automatically satisfied). 
Equations fig)). (fT5)l and ([To]), then lead to 
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FIG. 5: Critical temperature T,f for the F/S/F structure of 
Figure l.a in the P configuration (top panels) and difference 
AT C = T^ p — T c p between the critical temperatures in the P 
and AP configurations (bottom panels), as a function of the 
thickness cLf of the F layers. The left panels show the effect of 
a finite 7^ and the right panels the effect of a finite 7^. The 
four panels show the case 7^ = 7^ = with full red lines, 
for comparison. The other parameters used in the Figure 
are 7 = 1, j T = 0.12 and d s /£,s = 1- The SDIPS modifies 
T% and AT C in a quantitative or qualitative way, depending 
on the case considered. The type of effects produced by the 
SDIPS on To and AT C strongly depend on the signs of j F 
and 7^. 



A c 



with 



b AP = 



u n \ + 2&BCsb c n J s 1 
AC, 



8Re[£° n J+4 


C° 


2 6s 


8 + 6Re[£ Q na }S s + 


c° 

1 


,<7 


°s 



and 



7t(S^;2 + i 7 ^crsgn(w„)) 



7 7 t + Bn,a + crsgn(cj„) 
These Eqs. together with ([3]) lead to 



(22) 

(23) 
(24) 

n$0-sgn(w„) (25) 



'rpBCS 

log ( ~^c- 1 = Re 




b c T 



BCS 



2 ^exp(r) r c c 



(26) 



where T denotes Euler's constant. The resulting Tjf is 
independent from the values of n and a used in Eq. (|26[) . 
Note that in the case 7^ = j F = and <5s — * 0, this equa- 
tion is in agreement with Eqs. (18) and (19) of Ref. [35L 

Performing a numerical resolution of Eq. (|26p together 
with (|2"3")) and (|24[) . one obtains the results of Fig. 
which shows the critical temperature T F of the structure 
in the P configuration (top panels) and the difference 
AT C = T AP — T F (bottom panels) as a function of dp, 
for different interface parameters^. Here, for simplic- 
ity, I consider cases where Tjf does not show a strongly 
reentrant behavior, i.e. a cancellation of XjF in a certain 
interval of dp. Such a behavior can happen for instance 



but 



for larger values of 7 t (see e.g. Ref. |35|) or 7^ 
this corresponds to a minority^ of the experimental ob- 
servations made so far. The four panels of Fig. O show 
the case 7^ = "f F — with full red lines, for comparison 
with the other cases, where the SDIPS is finite. I first 
comment the results obtained for the T F '(dp) curves (top 
panels of Fig. [5]). In some cases, the SDIPS can modify 
quantitatively the T F (dp) curves, for instance by am- 
plifying the dip expected in T c p (dp) in the absence of a 
SDIPS (see e.g. dotted curve in the upper left panel, cor- 
responding to 7^ > 0). In some other cases the SDIPS 
can modify qualitatively the T F {dp) curves, for instance 
by transforming the minimum expected in T F (dp) into a 
maximum (see e.g. curve with circles in the upper right 
panel, corresponding to 7^ < 0). I now comment the 
results obtained for the AT c (dp) curves (bottom pan- 
els of Fig. [S]). In the range of parameters studied in 
this work, one always finds AT C > because the effects 
of the two F layers on S are partially compensated in 
the AP case. In the absence of a SDIPS and for a low 
7t, the AT c (dp) curve presents a maximum at a finite 
value of dp (sec red full lines in bottom panels). In some 
cases, the SDIPS can modify quantitatively the AT c (dp) 
curves, for instance by increasing the value of this max- 
imum (see dotted curve in the bottom left panel, cor- 
responding to 7^ > 0). In other cases, the SDIPS can 
modify qualitatively the AT c (dp) curves, for instance by 
turning the maximum expected in AT c (dp) into a mini- 
mum (see dashed curve in the bottom left panel, corre- 
sponding to 7^ < 0), or by increasing the complexity of 
the variations of AT C with dp (see curve with circles in 
the bottom right panel, corresponding to j F < 0), or by 
transforming AT c (dp) into a monotonically decreasing 
curve with a reduced amplitude (see dot-dashed curve in 
bottom right panel, corresponding to j F > 0). Remark- 
ably, the type of effects produced by the SDIPS on the 
T^(dp) and AT c (dp) curves depend on the sign of j F 
and 7?. Critical temperature measurements can thus be 
an interesting way to determine the values of 7^ and 7?. 
For simplicity, I have shown here results for 7^ ^ and 
7^ ^ separately. Nevertheless, the behavior predicted 
for 7^ ^ and 7^ ^ simultaneously remain highly in- 
formative on the values of 7^ and 7^. Note that if j F 
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and 7^ are increased compared to the values used in Fig. 
El the T^((1f) curve gets some cancellation points like 
in Ref. [33, but the behaviors of Tjf and AT C remain, in 
many cases, qualitatively dependent on the signs of 7^ 
and 7^ (not shown). 

Before concluding, I note that with the approximations 
used in the previous section (III.B), the electronic corre- 
lations inside S were affected by the presence of the F 
electrodes through the parameter 7^ only. Consequently, 
the self-consistent gap A c of the S layer was independent 
from cIf, and it was furthermore identical for the P and 
AP configurations for 7^ = 0. In the present section, I 
have not neglected the dependence of and T C AP on dp 
and 7^ because I have considered parameters for which 
hypothesis 3 is not acceptable anymore. In particular, I 
have used lower values for dp21. 



IV. CONCLUSION 



This article shows that the Spin-Dependence of In- 
terfacial Phase Shifts (SDIPS) can have a large variety 
of signatures in diffusive superconducting/ferromagnetic 
(S/F) heterostructures. Ref. [29| had already predicted 
that the SDIPS produces a phase shifting of the oscilla- 
tions of the superconducting correlations with the thick- 
ness of F layers. This article shows that this is not the 
only consequence of the SDIPS in S/F circuits. In partic- 
ular, the SDIPS can produce an effective magnetic field 
in a diffusive S layer contacted to a diffusive F layer. 
This effective field can be seen e.g. through the DOS 
of the diffusive F layer, with a visibility which oscillates 
with the thickness of F. The SDIPS can also modify 
significantly the variations of the critical temperature of 
a S/F bilayer or a F/ S/F trilayer with the thickness of 
F, either in a quantitative or in a qualitative way, de- 
pending on the regime of parameters considered and in 
particular the sign of the conductances Gf and G% used 
to account for the SDIPS of the S/F interfaces. In the 
case of a F/ S/F spin valve, this last result also holds for 
the thickness-dependence of the difference between the 
critical temperatures in the parallel and antiparallel lead 
configurations. These effects should help to determine 
the parameters G§ and of diffusive S/F interfaces. 
The calculations shown in this paper are also appropri- 
ate to the case of thin diffusive S layers contacted to 
ferromagnetic insulators. 

I thank P. SanGiorgio for showing me his experimental 
data prior to publication, which stimulated part III.B of 
this work. I acknowledge discussions with T. Kontos and 
W. Belzig. This work was supported by grants from the 
Swiss National Science Foundation and Region Ile-de- 
France. 



V. APPENDIX A: PARAMETERS OF A S/F 
INTERFACE FROM A DIRAC BARRIER MODEL 



= 0.8 eV, Ej , = 1 2 eV, P = a = 0.06 




FIG. 6: Conductances Gt (full line), |G^| (dashed lines), and 
\G^\ (dash-dotted lines) reduced by the number of channels 
n, as a function of the spin-averaged barrier strength Z of a 
S/F interface modeled with a Dirac barrier (see text). The 
curves were obtained with typical Fermi energies Ef = 0.8 eV 
and E* = 12 eV in F and S respectively, a spin-polarization 
P — 0.06 of the density of states in F, and a spin asymmetry 
a = 0.06 for the barrier. Note that for the set of parameters 
used this Figure, one finds G^ < and G^ < 0. However, for 
other parameters (e.g. by using a < 0), one can reverse the 
signs of Gj and G^, or obtain opposite signs for Gj and G^ 
(not shown). 



The exact values of the conductances Gt, G^ and G^ 
of a S/F interface depend on the details of this inter- 
face and on the microscopic structure of the contacted 
materials. Nevertheless, it is already interesting to study 
a simplified Dirac barrier model which shows that the 
parameters regime assumed in this paper is, in princi- 
ple, possible. Neglecting the transverse part of the elec- 

(fc c</ ip\ k 

1/2 

Here, k 
and 



trons motion, one finds Tn a 



and tn\P 



2(kgk F ) 



S(F) 

with D a 



V F(S) 

k 



-iZ a )/D° 
k% + iZ° 



/D° w^ - -as-rap 

s ,p>. is the Fermi electronic wavevector in S(F) 
Z a is the strength of the Dirac barrier for elec- 



hk a s 



{2m e Ef) 1 l 2 , with 



trons with spin a. I use 
m e the free electron mass and Ef the Fermi level in 
S. For F, I use an s-band Stoner model, in which 
hk F = (2m e Ef[l ± 2P/(1 + P 2 )]) 1/2 , with P the spin- 
polarization of the density of states in F and Ej the 
Fermi level in F . I assume that the barrier can have a 
spin dependence a = {Z* — Z^)/(Z^ + due to the 
ferromagnetic contact material used to form the inter- 
face. The results given by this approach are shown in 



G 



and 



reduced 



Fig. The conductances Gt, 

by the number of channels n, are shown as a function of 
the spin-averaged barrier strength Z = (Z* + Z' l )/2. The 
three types of conductances go to zero when Z goes to 
infinity with a constantly. One can see that in the limit 
T n < 1 (i.e. here G T h/ne 2 < 1) and P <C 1 in which the 
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boundary conditions (|5l6p have been derived, G 1 ^ 

be significantly stronger than Gp, as assumed in Figs. [5] 
and [3] Note that for the set of parameters used in Fig. 
IH one finds G F < and G^ < 0. However, for other pa- 
rameters (e.g. by using a < 0), one can reverse the signs 
of G F and G^, or obtain opposite signs for G F and G^ 
(not shown). This suggests that there is no fundamental 
constraint on the signs of G F and Gf in the general case. 



VI. APPENDIX B: EFFECTIVE FIELD EFFECT 
IN A F LAYER 

This appendix reconsiders the case of the S/F bilayer 
of Fig.l.b, in the limit g?f *C £f where 9 n ,a(x) can be 
approximated with a quadratic form in the F layer. From 
Eqs. © and ©, one finds, for x 6 [ds/2, d s /2 + d F ], 



and 



C = 



Isin^j-Acos^) 
E s 



(32) 



The notations used in the above equations are the same 
as in Section III. In particular, the Thouless energy and 
the normal state conductance of the S layer with thick- 
ness ds/2 are denoted Gs — 2as A/ds and E§, H — 
4HDs/d s , and one has r] F '^ — ersgn(w„). In the limit 

6 F la < 1, Eq. JSl]) is in agreement with Eqs. (TF|) 
and (fl~8"|) . The analytic continuation of Eq. (|3"Tj) shows 
that the 5* layer is subject to the effective field H e ff = 
E^ H G^/Gs, in agreement with Eq. (f20|) . The analytic 
continuation of Eq. (f2"5)) shows that the F layer is subject 
to an analogue effective field defined by 



ds — 2x f ds — 2x 



dp 



2d F 



(27) 



F W 



(33) 



with 



: arctan 



7x sin I 



7 T cos(^ ) 



2£ F n„ 



■ F P,2 
- 1% Vn.a 



^nSin« CT ) 



jpF 
rj TH 



9. 



iEe<cT)n'3 and e L 



el a {x = d s /2). Here, 
Gf = opA/dp and E F H = HDp/dp denote the con- 
ductance and Thouless energy of the F layer (note that 
dp <C £p E ex <C E F H ). From completeness, I also 
give the analogous equations for x G [0, ds/2]. Assuming 
that 6 nM (x) can be approximated with a quadratic form 
in the S layer, one finds, from Eqs. ([!]) and (0, 



9 s - 



2{2x - d s ) / 2x-d s 
ds \ d s 



(30) 



with 



77t sm( 



, a" ) 



arctan 



4gs _A 
rfs E s - 



77T COS( 



9 F 
n.er 



) + 



S P.2 



It is interesting to compare Eqs. (|28 
that in the regime 



9 f 

n.a 



and Q. First, note 
can be recovered 



< 1, Eq 

from the low dp-limit of Eq. ((9|) . The interest of Eq 
is that it is valid beyond the regime of a weak supercon- 
ducting proximity effect studied in part III. In particular, 

(28) Eq. (f28|) indicates that, in principle, the field E F H can 
occur in a thin F layer for any value of the tunneling 
conductance Gp of the S/F interface. Interestingly, in 

(29) the case of a thick F layer dp > £ F /2, from Eq. ©, the 
contribution of G F to the pairing angle F a (x) cannot 
be put anymore under the form of an effective exchange 
field, due to the non-linearity of the term. Thus, 
strictly speaking, the concept of a SDIPS-induced effec- 
tive field is valid only in the limit of thin metallic layers. 
However, it might be possible, in principle, to observe 
reminiscences of the double-gap structure appearing in 
N F (e) in the regime of intermediate layer thicknesses 
dp ~ £f at least, due to the continuity of the equations. 
Then, one can wonder why this effect does not appear 
in section III. This is due to the regime of parameters 
chosen: section III assumes E ex 3> A^cs like in most 
experiments, and it furthermore focuses on the typical 
energy range probed in superconducting proximity effect 
measurements, i.e. \s\ < 2Abcs- I n such a regime, one 
can neglect the term \u) n \ compared to iE ex asgn(u! n ) in 
Eqs. and (|28[) . so that Fl F ^ does not emerge in the 
model. Accordingly, signatures of the SDIPS-induced ef- 

(31) fective field H F ^ are not likely to be observed in practice. 
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